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Abstract 

The Fateev model is somewhat special among two-dimensional quantum 
field theories. For different values of the parameters, it can be reduced 
to a variety of integrable systems. An incomplete list of the reductions 
includes 0(3) and 0(4) non-linear sigma models and their continuous 
deformations (2D and 3D sausages, anisotropic principal chiral field), 
the Bukhvostov-Lipatov model, the N = 2 supersymmetric sine-Gordon 
model, as well as the integrable perturbed SU 2 (n) <E> SU 2 (p — 2) / SU 2 (n + 
p — 2) coset CFT. The model possesses a mysterious symmetry structure 

of the exceptional quantum superalgebras U q (D(2\l;ot)). 
In this work, we propose the ODE/IM correspondence between the Fa- 
teev model and a certain generalization of the classical problem of con- 
stant mean curvature embedding of a thrice-punctured sphere in AdS%. 



March 2013 



1 Introduction 



Broadly speaking the ODE/IM correspondence is a link between the theory of Integrable Models 
in two dimensions and the spectral analysis of Ordinary Differential Equations. The approach 
was originated by Dorey and Tateo in Ref. [I] from the observation that the thermodynamic 
Bethe ansatz equations for certain 2D CFT minimal models proposed in [2] coincide with an 
exact version of the Bohr-Sommerfeld quantization condition for ID anharmonic oscillator - a 
remarkable result due to Voros [3H5]. Very shortly the initial observation was generalized and 
proved in Ref. [6]. Later the ODE/IM correspondence was established for a large variety of mod- 
els (for review see [7]). However, during the next decade, all attempts to incorporate massive 
integrable QFT in the ODE/IM correspondence have failed. Since the work of Gaiotto, Moore 
and Neitzke [S], thermodynamic Bethe ansatz equations have emerged in different contexts 
of SUSY gauge theories, algebraic and differential geometry [9HT2]. Inspired by this progress 
Zamolodchikov and the author have established the ODE/IM correspondence for the quantum 
sine/sinh- Gordon model - the model which always served as a basis for the development of the 
2D integrable QFT [T3]. Recently the result of Ref. [T5] was extended to the Toda type QFT 
model [33]. It has become clear that the ODE/IM correspondence is one of the most impor- 
tant aspects of integrability in 2D QFT. Among open problems within the approach is how to 
incorporate the class of integrable non linear sigma models into the ODE/IM correspondence. 
In this work, we propose an example which aims to step in this direction. 
We begin with the following well known fact (see e.g. Ref. [T5]); 

Let £ 9>n be a compact Riemann surface with n marked points ( "punctures" ) and 

be positive numbers such that 2x(S 9 ) + X]r=i( a « — ^) =0. Then there exists a flat metric 

on E 5 n with conical singularities of angle irai at the i th puncture. The metric is unique up to 

homothety. 

In the case of a two-sphere with three punctures the theorem's condition reads as 

ai + a 2 + a 3 = 2 , (1.1) 

whereas its statement is somewhat trivial. Indeed introduce a complex coordinate z and define 
a holomorphic differential p(z) (dz) 2 on the universal cover of £0,3 by means of the following 
assignment for p(z): 



[z - z A ) 2 ~ a ^{z - z 2 ) 2 - a2 {z - Z 3 ) 2 ~ a * 

Then the flat metric reads explicitly 

(ds) 2 = y / p{z)p{z) dzdz . (1.3) 

Here p stands for the homothety parameter and Z{ labels the punctures. 

Consider now the problem of constant mean curvature embedding of £0,3 into AdS 3 . In 
this case, the Gauss-Peterson-Codazzi equation can be brought to the form of the modified 
Sinh-Gordon (MShG) equation 

d z d- z 7] - e 2 " + p(z)p(z) e~ 2v = , (1.4) 
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where the field rj defines the induced metric [TO lfTGHTS] 1 ! 

4 e 2r? 

(ds)Lc = ; ~ m i , w . (ds)o , (1-5) 
1 + H y/p(z)p{z) 

and H = const stands for the mean curvature. A suitable solution should be real and smooth 
as z 7^ Zi, and, if we want to preserve the amount of the Gaussian curvature localized at the 
punctures, it should satisfy the conditions 

n — \ log (p(z)p(z) ) = O(l) at z — > Zj (i— 1,2,3) and oo . (1.6) 

It turns out that the problem (11.41) . ( jl~6i) admits an important generalization. One can 
consider the MShG equation in the flat background of the punctured sphere E ,3 subject of 
more general asymptotic conditions 

rj = —2 log \z\ + 0(1) at z — > oo 

n = Irrii log \z — Zi\ + 0(1) at z — > Zi . (1.7) 

If 



< ^ < 2 (1.8) 

and 

-\<m*<-\ (2 -a,) , (1.9) 

then the solution of the generalized problem exists and is unique. We are not going to prove 
this statement here. Instead, we will accept it and argue that the solution of (11.41) . (jl.7p is 
related to a certain 2D integrable QFT model, where all the parameters aj, p and m,^ as well 
as the restrictions ( JHHJ) and (ll.9p . admit simple physical interpretations. 

To describe the relation between classical and quantum systems we shall need an important 
property of the problem (11.41) . (jl.7p . The essential point in the formal theory of the equation 
(ll.4p is an existence of an infinite hierarchy of one-forms, which are closed by virtue of the 
MShG equation only, 

{t0 2 n}n=l ■ diO 2n = 0. (1.10) 

They are usually normalized by the condition 

002n = ( (VPW )^ 2n (9 zV ) 2n + . . . ) dz+ (...) dz , (1.11) 

where dots in the first bracket involves terms with higher derivatives of d z r) and/or p(z). The 
one-forms are not single-valued on the punctured sphere and should be considered on the 
universal cover of S .3 caused by the presence of the multivalued function *Jp(z). However, the 
contour C depicted in FigJT], loops around each puncture and yjp(z) acquires the same value 
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Figure 1: The integration contour C for the conserved charges q 



2n-l- 



after the analytic continuation along C. Therefore the restriction of bj<m to C are single-valued 
and the integrals 



(\2n-l 



'1.121 



c 



are not sensitive to continuous deformations of the contour. A branch of multivalued function 
\/p{z) can be chosen in such a way that the conserved charges are real numbers, 



(\2n-l — l2n-l 



;i.i3) 



We now turn to the QFT model of our interest. It was introduced by Fateev in Ref. [19 
and governed by the following Lagrangian in the 1 + 1 Minkowski space 



c = ^ £(cw-(W) 



i=l 



2/x (e 1 



cos(«iV9i + a 2 y?2) + e 



COS 



(ai<£>i - a 2 <^2) ) • 



Here ctj are coupling constants subject to a single constraint 



;i.i4) 



2 2 2 



(1.15) 



In this paper, we shall focus on the case where 

a\ > , «2 > 



> . 



1.16) 



1 The derivation of (|1 .4[) is nearly identical to the ones from textbook examples of embedding into M 3 , § 3 or 
H 3 (see, e.g. Ref. [H]). However in that cases Eqs. ([1.4[ ) is substituted by d z dzT) + e 2v — p(z)p(z) e~ 2v = 0. 
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The parameter /i in the Lagrangian sets the mass scale, /i ~ [mass]. We shall consider the 
theory in finite-size geometry, with the spatial coordinate x in </?j = ipi(x,t) compactified on a 
circle of circumference R, with the periodic boundary conditions 

<Pi(x + R,t) = (pi(x,t) . (1-17) 

Due to the periodicity of the potential term in fll . 14[) in ifi, the space of states % splits on the 
orthogonal subspaces %k\MM characterized by the three "quasimomentums" kf. 

ipi^ipi + 2n/ai : | ^ kl MM ) -> ^ I ^fci.te.fcs > ■ ( L18 ) 
The QFT (I1.14p is integrable, in particular it has infinite set of commuting local integrals 

Ll< 



of motion I^-i) ^Li-i; 2n = 2, 4, 6, . . . being the Lorentz spins of the associated local densities 



^ = r ^ [ s (^^a) 2fe +...], (i.i9) 

^° i+j+k=n 

where d± = \{d x =F <9 4 ) and . . . stand for the terms involving higher derivatives of <fi, as well as 
the terms proportional to powers of /i. The constant was found in Ref. [20] 

Mn) = _±_ (2^(1 - 2n)) n _ t (2al (1 - 2n)) w _ . (2a§ (1 - 2n)) w _ fc 

i\ j\ k\ (2n-lf (4af) 1 " J (4a|) 1 "J (4a|) 1 " fc ' 1 j 



where (x) n stands for the Pochhammer symbol. Notice that the displayed terms in f ll . 19j) with 
the given C^j. set the normalization of I^Li unambiguously. 
Of primary interest are the /c-vacuum eigenvalues 

hn-i = I R) = iLUiki} I R) , (1.21) 

especially the /c-vacuum energy 

E = 2 h . (1.22) 

In the large- R limit all vacuum eigenvalues Im-i vanish except I\. The vacuum energy is 
composed of an extensive part proportional to the length of the system, 

£ = #£0 + 0(1) at R^oo. (1.23) 

One of Fateev's impressive results concerning theory (11.14)) is an elegant analytical expression 
for the specific bulk energy 
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The main observation of this work is that the vacuum eigenvalues can be expressed in terms 
of the classical conserved charges ( 11. 12ft : 



/i" 1 (h-\R£ ) = di qi (1.25) 

/i 1 " 2 " hn-l = dn <\2n-l (n = 2, 3, . . .) . 

Here d n are constants, independent of fcj and R. With the normalization conditions for fyn-i 
and described above, d n reads explicitly as 

(_-\\n-l 3 

dn = (2n) 2n ~ 1 Y[r{2(2n-l)a') . (1.26) 

i=l 

The parameters of the quantum and classical problems are identified as follows: 

a\ = ^ (i = 1,2,3) 

M = - (2m* + 1) , (1.27) 

whereas the relation between dimensionless parameter fiR and p is given by 

HR = 2p . (1.28) 

Notice that the quantum integrals of motion I^n-i are periodic functions of the quasimo- 
mentums, and one can assume that ki takes values within the first Brillouin zone 

-\<h<\- (1.29) 

Since I^n-i commute with the "charge conjugations" - the obvious symmetries </?j — > — </?j (z = 
1,2,3) of the Lagrangian, the eigenvalues (ll.2ip are even functions of ki. In a view of the 
identification (ll.27p . the inequality (I1.29P suggests a natural domain (jl.9p for the parameters 
mj. Although both the classical conserved charges and the eigenvalues I<m-\ are nonsingular at 
fcj = (rrii = — |), the small- 1 2 — Zj| asymptotic (jl.7p contains a subleading log log-term in this 
case (the so-called "parabolic point"). For this reason, the value = — | is excluded from the 
domain ( II. 9p . 

Let us recall now that the MShG equation constitutes a flatness condition for a certain 
s/(2)-valued connection A = A z dz + A s dz. The associated linear problem 

(d z -A g )* = 0, (d 2 -Az)* = (1.30) 

is of prime importance for solving (jl.4p . (jl.7p . In particular, certain monodromy coefficients 
for the linear system ( ll.30p can be regarded as generating functions for the set of conserved 
charges {q2 n _i}^ =1 . Thus Eqs.( ll.25"j) relate spectral characteristics of the linear problem ( ll.30p 
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with a vacuum spectrum of the local integral of motions from the Fateev model. This is an 
example of the ODE/IM correspondence^ 

Al.B. Zamolodchikov was probably the first who realized the main advantages of ODE/IM 
correspondence compare to the traditional approaches [21]; The ODE counterpart makes ex- 
plicit the analytic properties of relevant physical quantities considered as functions of the pa- 
rameters. Within the ODE/IM approach, the integrable model can be studied uniformly in 
different parameter regimes. In the case of massive QFT, this was explicitly demonstrated in 
Ref. [13] for the quantum sine- and sinh-Gordon models. In spite of formal similarity of the 
Lagrangians, the physical content of these models is very different from one another. However, 
the models can be treated uniformly within the ODE/IM approach. With only a few minor 
modifications one can jump from the trigonometric to the hyperbolic model. The story is re- 
peating itself - only little modifications a la sin/sinh- Gordon one are required to extend the 
above ODE/IM correspondence to the most interesting regime of QFT f)1.14p with 

al>0, aj>0, a 2 3 <0. (1.31) 

In this regime, the Fateev model provides a dual description of the "3D-sausage" - the integrable 
non-linear sigma model with three dimensional target space [19J. The sigma-model regime will 
be the subject of a separate publication. The main purpose of this work is to present evidences 
in support of the relations fll.25p - fll.28p in the regime (j!.16p . 

The paper is organized as follows. In Section[2]we give an accurate definition of the conserved 
charges (\2n-i an d recall their relation to the flat connection. In the next section we focus on the 
first conserved charge qi. Our analysis is based on the observation that as p — > the problem 
fll.4j) . fll.7l) is reduced to the problem of constructing a metric of constant intrinsic curvature 
on the punctured sphere. Using the results from the classical Liouville theory we derive the 
small-p expansion of qi. SectionH] deals with the p — > limit of the conserved charges q2n-i 
with n > 1. Finally, in Section[5]we extract the relations f ll.25p - p.28p from the comparison of 
the small-p expansions for (\2n-1 against predictions of the Conformal Perturbation Theory for 
the QFT 01.14p . We conclude with a few remarks on a rigorous proof of the proposed ODE/IM 
correspondence. 

In the theory of integrable models, an important role belongs to the concept of the Yang- 
Yang function(al) [22]. It is of prime interest to the ODE/IM correspondence [23]. In this 
work, we use the Yang- Yang function as an auxiliary tool only to derive the small p-expansion 
of the conserved charge qi. Because of its own significance, we present some facts concerning 
the Yang- Yang function for the Fateev model in the appendix. 

2 Conserved charges for MShG on the punctured sphere 

In this section we give detailed description of the conserved charges {<j2n-i}5JLi- 

2 In fact, the abbreviation "IM" can be understood as a shortened form of cither Integrable Models or 
Integrals of Motion. The author slants toward the second interpretation; Until now there is no any indication 
that the correspondence can be extended beyond the scope of relations between the spectral characteristics. In 
particular, the relation of the formal variables (z, z) and the space-time coordinates (x,t), as well as the role of 
the classical field rj itself in the quantum theory, remains completely mysterious. 
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2.1 From MShG to ShG 



The flat metric defines a complex structure on a punctured sphere £0,3 • Since the MShG equa- 
tion (jl.4p is covariant under coordinate diffeomorphisms, its form can be simplified by taking 
advantage of conformal transformations. First of all, we can use the Mobius transformation 

c= (»-«,)(»-«.) (21) 

{Z2 ~ Zi){z - Z 3 ) 

to convert the puncture's coordinates (z±, z 2 , z 3 ) to their standard values (0, 1, oo). This brings 
Eq. ffL4"j) to the form 



where 



d c dtf - e 2 ^ + P(QP(0 e" 2 ^ = , (2.2) 



V = V + \ log (| g) , P(0 = P 2 r~ 2 (1 - CP" 2 ■ (2.3) 



AC dC 

Then we consider the Schwarz-Christoffel mapping 



"(0 = / dC x/nO , (2-4) 



which transforms the upper half plane Qm(() > to the triangle (wi, w 2 , w^) in the complex 
w-plane (see Fig. [2]). The lower half plane ^sm(Q < can be mapped into the triangle reflected 
w.r.t. the straight line (w\, 102). If there exists a unique solution, it should respect all symmetries 
of the problem under consideration. In particular it should be invariant under the reflection 

v(CX)=v(CO • (2.5) 

Let D be a four-sided polygon with the vertexes located at = {wi, w 2 , w 3 , w 3 } whose 
adjacent sides are glued together as it is shown in Fig. [2] to form a topological two-sphere. The 
overall scale of the polygon is governed by the homothety parameter p in Eq. fll.2p . Without 
loss of generality we shall assume that p > 0, then 

In this way we reformulate the original problem as a problem of solving of the ShG equation 

dvdnrj - e 2 ^ + e" 2 ^ = (2.7) 
in the open domain D with the boundary 

<9B = U 3 i=1 { Wi = iBi} . (2.8) 
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w 3 



Figure 2: Triangle (w 1 ,W2,w 3 ) is a w-image of the upper half plane Sm(() > under the 
Schwarz-Christoffel mapping (12. 4p . The point W3 is a reflection of W3 w.r.t. the straight line 
(wi, w 2 ). The domain D is obtained from the 4-polygon (w±, w 3 , w 2 , W3) by the identification of 
the sides [101,103], [102,^3] and [wx, [u> 2 ,^3], respectively. 

The function 

fj = fj - \ log(PP) = 77 - \ log(pp) (2.9) 

should be real and smooth for all z inside D and possesses the following asymptotic behavior 
near the boundary, 

fj = 2k log \w - Wi\ + 0(1) at w-tWi, (2-10) 

with 

^ = 27^1 _ 1 (i = 1|2>3)< {2U) 
2.2 Definition of the conserved charges 

The ShG equation (12.71) possesses an infinite set of continuity equations in the form 

du>F2n = 9 w G2n-2 ■ (2-12) 

The functions (F 2n , G 2n -2) are conventional tensor densities. They can be described as follows. 
Let 

u= (d w f,) 2 - d 2 w f, , v = (d w f,) 2 + dlfj . (2.13) 

Then 

F 2n = U n [u] , G 2 n- 2 = e~ 2f > U n _x[v} - 2 8 n ,x , (2.14) 



8 



where C/ n [tt] are homogeneous (grade(w) = 2, grade(<9) = 1, grade([/ n ) = In) differential 
polynomials in u of degree n (known as the Gel'fand-Dikii polynomials [25]). 



Here 



UJu 



A 



2 n n\ 



(2n- 1)!! 



A n • 1 . 



1 ^2 



4 d z + u - \ d- 1 u' 



(2.15) 



(2.16) 



and prime stands for the derivative. Thus, 



U 


u] 


= 1 , 






u] 


= u , 




u 2 


u] 


= u 2 


-\u" , 


U 3 


u] 


= u 3 


- \ (u'Y 


u n 


u] 


= u n 


+ ... , 



(2.17) 



u 



u 



where the last line shows overall normalization of the polynomials. The continuity equations 
imply that 



(2.18) 



are closed one-forms. 

Let us consider the integral 



(\2n-l 



(2.19) 



where C w stands for some contour in the open set D. Generally speaking, the integral depends 
essentially on the choice of integration contour. However, for a non-contractible loop such that 
the local densities come to themself when they are translated along C w , the integral is not 
sensitive to continuous deformations of C w . In this case c\2n-i can be treated as a conserved 
charge associated with the closed one form u 2n - 

In order to choose a suitable integration contour in (I2.19p . it is useful to make the change 
of variables and return to the original coordinates z and z. Using the relation 



the integrand in (12 .191) can be rewritten in terms of 

u = (d z n) 2 - d 2 z n , 
which is a single-valued field on £0,3. F° r example, 

App" — hp' 2 



u = p 



u + 



lQp 2 



(2.20) 



(2.21) 



(2.22) 
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and therefore, for n = 1 the integral (12.191) can written as 



qi 



c i 



dz ( ^— + — p a Upp" - hp' 2 ) 
\y/p 16 



dz y/p 2 (^pp e~ 2v - 1) 



(2.23) 



Here C is an image of C w under the inverse conformal transformation w — w(z). Of course, 
it is straightforward to perform the change of variables for any given n. We do not present 
explicit formulas for n > 1, just notice that 

F 2n [u]=p~ n F 2n [u} + ... , (2.24) 

where the omitted terms contain derivatives of u(z). 

We now note that ^/p(z) acquires the same value after the analytic continuation along the 
contour depicted in FigJlJ Therefore the integral (I2.23P does not depends on a base point (an 
initial point of the integration path), as well as the precise shape of the contour. The image 
of this contour under the Mobius transformation (12.11) is the Pochhammer loop Cq shown in 
Fig. El Recall that the Pochhammer loop, considered as an element of the fundamental group 
7Ti(Eo,3), is given by the commutator of the loops 70 and 71 which wind anticlockwise around 
the points ( = and ( = 1, respectively: 



a 



7o 7i 7o 1 7i 1 



(2.25) 



Also, the complex conjugated contour is homotopically equivalent to the C^, 

Q = 7i~o 1 C C 7io w ith 710 = 7i 7o • 



(2.26) 



This way we can define the conserved charges by Eq. ( 12.191) where C w stands for the w-image 
(up to homotopy) of C$ (see Fig.[3j). One can introduce another set of conserved charges 



where 



(\2n-l 



3 i7r(n-i)(ai+a 2 ) 



( a 



dw F 2n + dw G 2n -2 



(2.27) 



(2.28) 



and the local densities (F 2n , G 2n - 2 ) are differential polynomials in 17 of the degree 2n which are 
obtained from (F 2n , G 2n _ 2 ) ( 12. 14|) by a formal substitution d w — > d^. The integral is taken over 
the contour C w which is complex conjugate to C w . 

Some explanation is in order here. Our definitions of the conserved charges suffers from 
an intrinsic phase ambiguity of the form e 17rn ( Ara i+ Ma 2) ^jy^ M G Z) inherited from the phase 
ambiguity of y / p(z). In what follows we assume the condition 



> 



as 



X\ < z < X 2 < X3 



(2.29) 
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w 



w, 




Figure 3: The Pochhammer loop CV with the base point (q. Its w-image is homotopically 

f =* -± T» 

equivalent to a union of the oriented line segments C w = wqO U be U cd U da U au>o- The points 
{a, 6, c, d} in the £-plane and their images are denoted by the same symbols. 



and define \Jp(z) outside of this real domain by means of the analytic continuation. In terms 
of the coordinate ( this convention implies that the principal branch of \VP(Q in Eq.f l2.4p is 
chosen to satisfy the condition 



y/P(0 > as < C < 1 • 



(2.30) 



For practical calculations, it is convenient to rewrite the defining relations (12.191) . (I2.28P in terms 
of the (C, C)-coordinates. Then the integrals should be taken along the Pochhammer loops 
and Cq- Bearing in mind condition (12.301) . it makes sense to chose the base point Co G C*c 
within the segment [0, 1] as shown in Fig.[3j The values of the integrands are determined 
unambiguously through the analytic continuation along the integration path. This convention 
resolves the issue of phase ambiguity in the definitions of (\2n-1 and c\2n-i- It also implies that 
the local densities (-F 2n , G2 n -2) are complex conjugates of (F 2n , G2n-2) 



2n 



2n ) 



G2n-2 — G* 2n -2 ■ 



(2.31] 



Now one can understand the role of the phase factor e 17r (™-i)( a i+ a 2) j n ^ e definitions of the 
conserved charges. Indeed, one easily verifies that Eqs. (l2.26p and (I2.3ip yield the complex 
conjugation condition 



l2n-l — l2n- 



(2.32) 



In the case under consideration the reflection symmetry w <->• w of the ShG equation is not 
broken (i.e. f)(w,w) = f)(w,w)), and therefore 



t?2n-l — l2n-l • 

For this reason, we can focus on the integrals (\2n-1 only. 



(2.33) 
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2.3 A generating function for the conserved charges 

The MShG equation constitutes a flatness condition for the sl(2) connectiorJl 

A z = -\d z r]a z + e e [a + e" 1 + a~ p{z)e- ri ] (2.34) 
A- z = \ d- zV a 3 + e~ e [a~ + a + p(z) e"" ] , 

with the spectral parameter 9. The connection A = A z dz + A z dz is not single-valued on the 
punctured sphere. However, it does return to the original branch after the continuation along 
the non-contractible loop C depicted in Fig.[TJ Therefore the Wilson loop 



W(9) = Tr 



Pexp (j> A ) 



(2.35) 



does not depend on the precise shape of the cycle used. 

Apparently, W(9) is an entire function of 9. Since the shift of the argument 9 — > 9 + in can 
be compensated by the similarity transformation A — > a 3 Aa 3 , the Wilson loop is a periodic 
function of period in 

W(9 + ivr) = W{9) . (2.36) 

It is easy to see that 

W(9) ^2cosh(q e e + o(l)) as Ue(9) -> +oo , \$tm(9)\ < | , (2.37) 



where 



lo = j dz^fpjz) 



q = <b dz^p(z) . (2.38) 



The integral does not depend Zi and we can set them to the standard values (0, 1, oo). Then 
the integral over the Pochhammer loop is performed using well known relation 

dC ^(1 - C)^ 1 = (1 - ^) (1 - e 2 ^) , (2.39) 



yielding the result 



qo = - 3 4?rW • (2.40) 



The subleading terms in the asymptotic expansion defined by the conserved charges. The 
textbook calculation yields the following asymptotic expansions |26j : 

, wm I -qoe e + Er=i^q 2 „-ie-( 2 - 1 ) e as Ue{9) +oo, |3m(0)| < f 

-qoe- e + E"=iCnq 2 n-ie^- 1 ) e as Ue{9) -+ -oo, \%m(9)\ < f ' 1 ' 1 



3 er- axe the usual Pauli matrices, i.e., a* = ( J , = (J J) , a" = (J J 
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here c n = ^7^. Thus W{9) can be regarded as the generating function for the the 
conserved charges. 

Notice that in the case of unbroken reflection symmetry the Wilson loop is an even function 
of the spectral parameter 

W{6) = W(-0) . (2.42) 

3 Small p-expansion of the first conserved charge 

Among the conserved charges, qi and cji are of special interest. Their linear combinations 

& = qi + qi , p = qi - qi . (3.1) 

can be understood as the energy and momentum of the Pochhammer string. In the case under 
consideration, the momentum is zero. Here we aim to explore the small-p expansion of qi. 

3.1 On-shell value of the action functional 

The ShG equation ( 12. 7p as well as the asymptotic conditions f)2.10p follow from the action 

- / d 2 w ( dwfjdujf) + 4 sm 2 (fj) ) + — / d£fj — log(e) 1 
. * Jul**) ~( ne J a ~( 

The first term here contains a two-fold integral over the domain 

D( rc g) = D\ U? =1 , B® = {weB:\w-Wi\<e} (3.3) 

depicted in Fig.HJ It is the "cutoff" version of the naive action functional. The additional 
terms involve the integrals over the boundary cffi)( reg ) (a union of the three closed contours C±, 
C2 and C3 UC3) and field-independent "counterterms" which provide an existence of the limit. 
We define the on-shell action A* as a value of the functional A[fj] calculated on a solution of 
the problem (l2T7jl . f l2~T0|) . 

It was already mentioned that the overall size of the polygon is defined by the homothety 
parameter p. As p — > +00 the dominating contributions to the on-shell action come from 
the boundary dB>. Near the vertexes, 17 can be expressed in terms of certain Painleve III 
transcendents. This observation allows one to determine (see, e.g., Appendix B from Ref. [23] ) 
the limiting value of the on-shell action 

Al = lim A* . (3.4) 

p— >+oo 

The result reads as follows 

3 

■A*oo = -J2 a > 0* log(2) + F + , (3.5) 



A\fj] = lim 



Qi I 2 



(3.2) 



'in this work we always use the shortcut notation for the 2D integration measure: d 2 w := ^ dw A dw. 
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c 3 

Figure 4: The integration domain B( reg ) for the regularized action (13.21) . 



where we use the notations 

F(x) = 



dt log 7(2) 



j(x) 



T(x) 



r(i -x) 

It makes sense to subtract the constant fl3.4[) from the action f)3.2p and define 

This function is the on-shell action normalized by the condition 

lim 2) = . 



(3.6) 



(3.7) 



(3.{ 



We shall refer to 2) = 2J(p|{/j, CKj}) as to the "Yang- Yang function" (for general explanation of 
this term, see Ref. [22]. In this work we closely follow the notations from Ref. |23j). 

3.2 Small-/? expansion of the Yang- Yang function 

To explore the small p-limit it is useful to return to the original complex coordinates (z, z) and 
the field rj. The action (13. 2 j) can be rewritten in the form (see Appendix A for details): 



A[ V ] 



Ao + lim 



7T )\z-*i\>t 
\z\<R 



d 2 z ( d z r]d- z r] + e 2 " + p{z)p{z) e" 2 " ) 



(3.9) 



+ 2 ^(m^i-m 2 log(ei)) +2^ + 2 \ogR 



Here 



27re, ; 



d£r) 



\Z-Zi\=£i 



2txR 



dlrj 



(3.10) 



\z\=R 
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whereas an explicit form of the field-independent constant Aq is given by the equation ( 1A.5I) 
from Appendix A. 

In the vicinity of punctures one has 

e 2 " oc \z - Zi\ Ami , p(z)p(z) e~ 2ri oc p 4 \z - Z .| 2 K- 2 ^- 2 ) . (3.11) 

Therefore, assuming that 

- \ < mi < -§ + \ en , (3.12) 

we may ignore the term p(z)p(z) e~ 2ri in the functional ( 13. 9 j) as p — > 0. In other words, the 
small-p limit is controlled by the classical Liouville theory on the sphere with three punctures, 
and the following limit does exist for any z ^ zc 

\imrj(z,z) = \ (p(z,z) . (3.13) 

Here the field tp is a solution of the Liouville equation 

d z d g ip = 2&, (3.14) 

subject of the boundary conditions 

<f(z,z) = Ami log \z — Zi\ + 0(1) at z — > Zi 

<f(z,z) = -4 log \z\ + 0(1) at z -> oo . (3.15) 
This is exactly the problem of constructing a metric of constant intrinsic curvature 

ds 2 cic = e v dzdz (3.16) 
on the sphere with three punctures. The total intrinsic curvature is given by 

r 3 

/ %]=47r (l + XX) . (3.17) 

Since m, obey inequality (I3.12p . then — 2ir < J S[)a K < 0. Therefore the solution exists and 
defines a metric of constant negative Gaussian curvature on £0,3. The corresponding on-shell 
value of the Liouville regularized action was calculated in Ref. [23]. Using this result it is 
straightforward to derive the small-p behavior of the Yang- Yang function ( 13. 7\i (see also the 
end of Appendix A for some technical details): 

3 3 
2) = MP) J2{l- ± P')+^- 2 P 2 U^) +0 ^- ( 3 - 18 ) 

i=l a% i=l 

Here 

Pi = mi + \ (3.19) 
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and the p- independent constant is expressed in terms of the function F ( 13. 6ft . 

2)0 = ^(3 -Vi -P2-Ps)+F(l - Pl -p2+P3)+F(l-p 2 -p3 + Pl ) (3.20) 

3 

+F{± -P3- Pl + P 2) - F(0) +52(aiF(%)-F{l- 2 Pl ) + a t (3a - |) 2 log(2) 



i=l 



Eq. ( 13.181) shows the nature of singular behavior at p — > 0. Notice that the term p(z)p(z) e~ 2v 
in the action ( 13. 9 p can be treated as a uniformly bounded perturbation oc p 4 , as far as the 
condition (13 . 1 2[) is fulfilled]^ Thus, the Yang- Yang function (13. 7p admits the small-p expansion 
of the form 

3 . 3 00 

?) = log(p) EU"^0 +?)o - 2p2 n^(|)+E^ ■ ( 3 - 21 ) 

i=l 1 i=l n=l 

3.3 Infinitesimal homothety of the flat metric of punctured sphere 

The variation of the action under infinitesimal perturbations of the moduli of the flat metric on 
£0,3 do not vanish on-shell. They can be expressed in terms of the stress-energy tensor. Here 
we focus on the infinitesimal homothety. 

For an infinitesimal dilation - = - = A<l, one has 



p 

Sp 



2 r 3 

lim - / d^B-Vfli lj , (3.22) 



where B stands for the trace of the stress-energy tensor of the ShG equation, 

6 = 4 sinh 2 (r/) . (3.23) 

The calculation outlined in Appendix B yields the formula which expresses the difference in 
the bracket in Eq. fl3.22|) in terms of the conserved charge qi = cji and the circumdiameter d of 
the triangle (u>i, u> 2 , ^3) in Fig.[2j 

p fr = ' (3 - 24) 



with 



i=\ 



Eq. (]3.8p implies that 5 is normalized by the condition 



lim $ = . (3.26) 

p— >+oo 



5 It seems likely that the formal proof of existence and uniqueness solution of the problem (|1.4j) . (| 1 . T[) can 
be obtained basing on this observation. 
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It is easy to see that 5 can be represented in a form 

3 3 



5=-i£(l-^)+Vn^(!)-^ / d^(*)e"\ (3.27) 

1=1 % 1=1 n ^^0,3 

and the small- p expansion (13.211) yield^H 

- / d 2 zp(z)p(z)e-^ = \2ny n p in . (3.28) 

The first coefficient in this series is simply expressed in terms of the Liouville field p (13.131) . 

Vx = — A I d 2 z p(z)p(z) e-*M . (3.29) 
vrp 4 J 

An explicit analytical expression for tp can be found in Ref. [21] (see Eqs.(4.1)-(4.5) therein). 



4 Small-p limit of Cfan-i f° r n > 1 

Eqs. fl3T25D -f l3T28D imply that 



It is instructive to verify this result using the original definition (12.23)) . 

The limiting value of the field rj is expressed in terms of the solution of the Liouville equation. 
Therefore as p — > the field u(z, z) turns to be a holomorphic component of the Liouville stress- 
energy tensor: 

\imu{z,z) = T L (z) with T L (z) = - (d z¥ >) 2 - \ d 2 ip . (4.2) 
In the case of the sphere with three punctures ( 13.151) the holomorphic fields Tl(z) has the form 

^) = -E((^ + ^-). where * = I-rf. (4.3) 



Since z = oo is a regular point, Tl(z) oc s as z — > oo. This imposes three linear relations on 
Q which determine them uniquely, 



<Sj + 5j - 4 - Sj 

Zj Z% Zfc Z{ 



{i,j,k) =perm(l,2,3) . (4.4) 



^Notice that, as it follows from the MShG equation, 



- / d 2 ze 2 " = - / d 2 zp(z)p(z)c-^ -lya t h 
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Combining Eqs. ( 12.231) and ( 14. 2ft . one gets 



qi 



dz 



c 



TAz) + 



1 8%p(z) 
4 p(z) 
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dzp(z) 
p{z) 



0(1) 



as p — > 



(4.5) 



where the contour is shown in Fig.[TJ Let Z\ < z% < z%, then the base point of C should be taken 
inside the interval [z± , z^\ , whereas the branch of the multivalued function yp{z) is determined 
by the condition (12.291) . With this convention, the integral does not depend on Zi and we can 
set them to be (0, 1, oo). Then using relation (12.391) . one easily reproduces Eq. (14.1j) . 

For n > 1, the leading small-p behavior of c^n-i can be obtained similarly. The calculations 
are elementary but rather cumbersome. For example, for n — 2, the result can be written in 
the form 



lim( p 3 q 3 ; 



p->0 



2tt 2 



r 3 



+ E^ 



1 

24 



i=l 



16 ~ 16 + 192 



4 



1 

24 



+ 



1 3 

240 ^ 

i=i 



where the numerical coefficients Ei, Eij and Hi are given by 

Ei = a,i (3aj - 2) (3a k - 2) 
2) 



E^j — 3 ct j (3(3^ 



#,: = 8 



9 (aiQ2 + «2fl3 + «3«i) + 15 aia2«3 



4p| 



(4.6) 



(4.7) 



In these equations (i, j, fc) represents any permutation of the numbers (1, 2, 3). It is much easy 
to establish the following general structure, 



lim( p 



i „ x _ o / 4p? 4p| 4p2 

H2n-1 J — -Kn > > 

ai a 2 a 3 



4„2 

where i? n stand for polynomials in the variables Xi = — of degree n, 



i? n (xi,x 2 ,x 3 ) = E ^i* 

i+j+k=n 



(n) + 



(4.9) 



The dots here represent the sum of monomials of degree lower than n. It is not difficult to 
calculate the highest coefficients for any n, 



R 



(n) 



-1 



in— 1 2n^-2 



n 



i! j! k\ (2n - l) 3 a}-*^^-* 



(4.10) 



18 



5 Identification of the parameters 



We are now in position to relate the parameters of the problem ( 11.41) . (jl.7p and the couplings 
of the Lagrangian (11.141) . As it was mentioned, the Fateev model has infinitely many local 
integrals of motion. The displayed terms in Eq. ( 11. 19ft fix the normalization of these operators. 

Let — I^n-i = ^2n-i be the vacuum eigenvalues of I^n-i aric ^ ^2n-i- I n ^ ne CFT limit, i.e. 
at fi = in the Lagrangian (j!.14p . these functions become polynomials in fcj of the degree n0 
and the normalization in (I1.19P is such that at n = we have 



'2n-l \n=0 



2?r\2n-i 



Cg2 (2« 1 fc 1 ) 2i {2a 2 k 2 f {2a 3 k 3 y 

i+j+k=n 



(5.1) 



The polynomials are known explicitly for n = 1, 2, whereas the constant is known for 
any n and given by Eq. (ll.20p . All these results were obtained in Ref. [20J. In a view of 
Eqs. (l4.ip . (I4.6p - (l4.10p . they are all in agreement with the proposal (11.251) . provided that the 
identification (I1.27P and the relation involving the normalization constant d n , 

d n 3 ( " ir " ^ limf^f-^l, (5.2) 

nti r( 2 (2n - 1) a?) ^oUnpJ v ; 

are imposed. 

To find the relation between p and the dimensionless parameter pR, let us consider the 
small- R expansion of the finite- volume energy (ll.22p . A brief inspection of the Lagrangian 
(I1.14P reveals that 

— = -f -jte n (»Rr. (5-3) 

7T O L — ' 

n=l 

Here the first term is dictated by the simple Gaussian model underlining the CFT limit with 
the effective central charge 

c cff = ]T(l-24c^ 2 ) . (5.4) 

The general structure of the short distance expansion follows from the fact that the potential 
term in the Lagrangian ( I1.14p is a uniformly bounded perturbation for any finite value of the 
dimensionless parameter pR. Therefore the Conformal Perturbation Theory can be applied 
literally and yields an expansion of the form (15. 3p with coefficients e n are expressed in terms 
of certain 2D Coulomb-type integrals. The large- R behavior of (15.31) is defined by the specific 
bulk energy (ll.23p . (11.241) . Eqs. ( I3.26]) -( I3.28|) strongly suggest the following identification 

d=- {E-RSo) • (5.5) 

7T 



' In this limit, (jl.ISp acquires continuous symmetries with respect to any shifts of the fields (pi, and the 
limiting values (|5.1[) are no longer periodic in ki. 
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This is, in fact, the first line in ( II .251) . provided that | pR coincides with p and the coefficients 
e n in expansion ( 15. 3ft are simply related to the coefficients y n from the small-p expansion of the 
Yang- Yang function ( !3.18p FI 

e n = 2 2 ~ 4n ny n . (5.6) 

Finally, Eq. (l5.2j) yields the formula ( 11. 26ft for the coefficient d n . 

The p — p relation fl 1 . 28[) implies that the bulk specific energy is simply expressed in terms 
of the area A of the punctured sphere S 3 calculated w.r.t. the flat metric (II. 3p (see Eq. (lA.6j) 
in Appendix A): 

R 2 S = -4A. (5.7) 

In order to explain the meaning of the lengths of the sides \wi — Wj\ (I2.6P in the Fateev model, 
let us recall some facts concerning the factorizable scattering theory associated with this QFT. 
All the details can be found in Appendix F in Ref. [27] . 

The spectrum consists of three quadruplets of fundamental particles 

Z% , e,e' = ±, z = l,2,3, (5.8) 



with the masses 



1=1 

and their bound states. (Here the relation = 4af is assumed to hold.) The Zamolodchikov- 
Faddeev commutation relations for the fundamental particles read 

£ 3 e 3 

Z^Z^) = ee"Y,[Sa k (0i ~ fc)]^ ^k^)^(^i) , (5.10) 

where (i,j,k) = perm(l,2,3) and 

S a (6)=i tanh(f + if ) S a (f + if ) . (5.11) 
Also S a (8) stands for the conventional 5-matrix in the quantum sine-Gordon theory [28] with 

(5.12) 



the renormalized coupling constant a, related to the Coleman coupling 0q [29] as follows 



J? 



2 



It is easy to see that the p — p relation implies that the dimensionless parameter | RM is 
the circumdiameter of the triangle Wiw 2 w 3 from Fig.|21 whereas jRMi are given by the lengths 
of the corresponding sides: 

RMi = 4 \wj — Wk\ , where (i,j,k) = perm(l,2,3) . (5.13) 



8 For n = 1, Eq. (j5.6p combined with p.29|) leads to a non-trivial prediction for the coefficient e\. It would 
be interesting to confirm this result within the Conformal Perturbation Theory. 
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6 Concluding remarks 



This work did not aim to achieve rigorous derivation of the ODE/IM correspondence. The 
goal was to propose a meaningful relation between a certain problem for the MShG equation 
on the one hand, and the Fateev model on the other. Forwarding can be performed along the 
following line. 

Having at hand Eqs.f l5.lOT) . it is not hard to guess the Non Linear Integral Equations (NLIE) 
[301 [31] associated with this factorizable scattering theory. In fact, the system of NLIE was 
already proposed by Fateev in his unpublished work [32]. Unfortunately, it still lacks the first 
principle derivation, say, from the latticed Nevertheless it would be important to confirm 
Fateev's NLIE from the ODE side. The part of the work was already done by Bazhanov. 

In order to explain the main result of the unpublished paper [31] , let us consider the auxiliary 
problem (ll.30p . (12.34]) . As is well known, this matrix system can be reduced to second order 
linear differential equations. One can write general solution of f l 1 . 3 j) as 

* = ( ei V ] = [ e ~ ' e ^ ( 9 * + d J^ $\ (6 l) 

\e"2 e~2 [d z + d z r]) ip I \ e"2 i/) I 

where if) and if) solve the equations 

[d 2 z - u(z, z) - e 29 p(z) ]if) = (6.2) 
[d 2 - u(z, z) - e- 2e p(z) ] ip = . 

Let us focus on the first equation in (16. 2p . This form is convenient for taking the limit p — > 0, 
when the field u(z, z) turns to be the holomorphic component of the Liouville stress-energy 
tensor Ti{z) (14. 3p . The function p(z) (ll.2p contains the small parameter p 2 as an overall factor, 
which can be absorbed by a suitable shift of the spectral parameter. Thus we define the new 
parameter 

k = p e e (6.3) 
and will keep it fixed as p — > 0. This yields the ODE 

[ -d 2 z + V (z) + V 1 (z)] 4 = 0, (6.4) 



with 



Vo(z) = -J2 



i=l 



Zi) 2 



VAz) - r 2 {Z3 ~ Z2r {Zl ~ Z3r2 {Z2 ~ Zl)a3 (6 51 



9 In the limiting case a§ = (the so called Bukhvostov-Lipatov model) the NLIE were derived from the 
coordinate Bethc Ansatz in Rcf. [551. 
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subject to the following constraints imposed on the parameters 

8i + 8j -8 k 8i + 8 k -8j . . . 

ai + a 2 + a 3 = 2, q = ^ 1 (i, j, k) = perm(l, 2, 3) . 

In the case k = the equation turns out to be Riemann's differential equation. By the simple 
change of variables and parameters, the ODE (16.41) can be reduced to the form used in Ref . [2U] . 
Some particular cases of this equation were studied in a series of works on integrable models of 
boundary interaction [351437] . 

Eq. ()6.4p was the starting point of the work [31]. Bazhanov derived a system of integral 
equations for certain connection coefficients of the ODE (16.41) . It occurs to be identical to 
Fateev's NLIE taken at the CFT limit. Needless to say that the limiting form of the NLIE 
differs from the general one by the source terms only. In all likelihood the original Bazhanov 
derivation can be adapted to the massive case. 
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A Derivation of Eq.(l318f) 



The purpose of this appendix is to rewrite the action functional (13. 2p in terms of the field rj 
and coordinates (z, z). We also outline the derivation of Eq. ()3.18p . 

Under the conformal map w — > z the image of the arc Ci enclosing the vertex Wi is an 
infinitesimal circle of radius e i: related to the cut-off parameter e ( 13. 3 p as 



2p 



Zjk 



2 Oj 
^ 2 



(a.i; 



Here (i,j,k) stands for any permutation of (1,2,3) and \zy \ = a/ (z; — Zj) (zi — Zj). One easily 
verifies the relation 



1 

TiC 



d£r) = aii]. 



c, 



y log 



Zjk 



with 



and then 



1 



cii(2 - a, ; 



log(ei) - di log(p) 



(A.2) 



d£r] 



Tie 



?-4 



i J\z-Zi\=ei 



2 - di (2 - a,j) 2 
d£r)-ai Z? log(e) = 2 (mi rji - m- log(ei) ) H ^ rji H ^— log(ej 



2mi(mi 



) log 


Zjk 




Zij^ik 



diliik + 1) log(p) + Cli If log 



(A.3) 



22 



The remaining necessary ingredient is 



lim 

e;^0 7T 



lim 

e e ->0 
R— >oo 



/ d 2 z d z f\d- z f] = ~\^{a 1 i- 4) log 

J\z-Zi\>ti ° i=1 

- [ d 2 z d z r)d 2 r) ~y^( ??i + 

7r / u-h\>h V 2 

|z|<H 1=1 



+ 



(2 - a,) 2 



(A.4) 



log(ei) ) + 2 r ]oo + 2\ogR 



where 



2ixR 



dirj . 



\z\=R 



Combining it with (1A.3[) . one arrives to Eq. fl3.9p . where the constant Aq is given by 



t=i 



+ E( i_(M )logW + Eoii?log (|)-^, 



t=i "* t=i 

and A is the of S 0j 3 w.r.t. the flat metric ( 11. 3ft . 



/. 3 

A = / d 2 2; ^fp{z)p(z) = np 2 J] 7 ( 

i=l 



The on-shell value, Al iouv , of the Liouville regularized action 



A 



Liouv 



lim 

£->0 
R— >oo L 



(A.5) 



(A.6) 



3 

/ d 2 2 ( <9 2 <^ + e^ ) + V ( m, ^ - 2 m 2 log(e) ) 

|z|<-R i=l 



where 



2 logi? + ^ c 



1 



(A.7) 



2?re 



V9j = * d£ If , V9 C 



|z— Zj|=e 



27T.R 



d£yj 



\z\=R 



was calculated in Ref. [23]. The authors found that the quantity 



1=1 



( 5- = -m»(m» + 1) ) 



(A.8) 



can be expressed in terms of the function F ( 13. 6 p ad_| 

g(cl) _ p^_ mi _ m2 _ m3 _ ^ _|_ p^_ mi _ m<1 _|_ m3 ~) _|_ p[^_ m<1 _ mg _|_ 

+F(-m 3 - mi + m 2 ) - F(0) - F(-2mi) - F(-2m 2 ) - F(-2m 3 ) . (A.9) 
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See Eqs.(3.21), (4.12) in Ref. [21] where one should set irjib 2 = 1 and rji is replaced by (- 
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Combining this result with ( 1A.5I) . one founds 

A* = S**> +J> H log (I) + £ (1 - ^i±D!) l 0g(p) - ^. a + o( P ) (A.10) 

i=l i=l 4 

as p — t- 0. Finally, using the formula (13. 5 j) for the constant .4.^, one obtains Eq. (13. 18[) . 



B Scalar potential for the stress-energy tensor 



Here we aim to obtain some useful identities involving qi and cji. 

The conserved charges qi and cji can be expressed in terms of the conventional stress-energy 
tensor associated with the ShG equation 



qi = e f (ax+oa) 



I (dwT + dwQ), q 1 = e ^ {ai+a2) I (dw T + dw 6) . (B.T 



Here T, T and stand for the non-vanishing components of the stress-energy tensor 

T=(d w fi)\ T=(<9^) 2 , 6 = 4 sinh 2 (7)) . 
The continuity equations 

d^T = d w Q , d w T = d^Q 
imply that the fields (IB.2p can be expressed in terms of a local potential function, 

T = <9?,$ , T = d 2 T ,^ , 6 = SUM 

and therefore, 



w ' 



qi = e 



f (ai+o 2 ) 



® (dw (9^ + dw du,) d w $ 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



The integration contour here can be chosen to be a union of oriented segments as shown in 
FigJHJ It is straightforward to express the integral as a sum of the discontinuities 



Ax 
A 2 



2i 
2i 



it; ^ 1 6 



i7ra2 i7ra2 

e 2 <9,„$L - e~&„$| 



as 6 G [wi, W3] 
as a G [u> 2 , w 3 ] . 



(B.6) 



Since qi does not actually depend on the location of the points a, c G [w) 2 , W3] and b, d G [u>i, W3], 
A x and A 2 remain constant along the sides [w 2 ,w 3 ] and [wi,w 3 ], respectively This yields the 
relation 



qi = sin (22a) Ax + sin (^f ) A 2 . 



(B.7) 
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Similarly one finds 



qi = sin (22a ) A{ + sin (^f) A* . (B.8) 

Of course, in the case under consideration Ai and A 2 are real constants. 

The formula flB.7j) allows one to simplify the w.h.s. of Eg. 03.22]) from the main body of 
the paper. Indeed, the on-shell value of the trace of the stress-energy tensor is given by the 
Laplacian of the scalar potential $ ( 1B.4I) . Hence, 



02) 



3 

lim - f d 2 w d w d^ - V at L 2 . (B.9) 



As it follows from the boundary conditions (I2.10p . 

$(w,w) = -2/ 2 log \w - Wi\ + 0(1) at w^Wi. (B.10) 

Combine this relation with the above observation that the discontinuities Ai and A2 (IB. 61) 
remain constant along the corresponding sides of the polygon (wi,w 3 , W2, W3), one finds 

<92J 1 

p — = -— ( Ai \w 3 - w x \ + A 2 \w 3 - wi\ ) . (B.ll) 
op 2tt 

Using (IB. 7j) . the w.h.s. of (IB. lip can be rewritten in terms of the conserved charge qi = qi 
and the circumdiameter d of the triangle (wx,w 2 ,w 3 ) (see Fig. [2] and Eq. (l2.6j) ). This yields 



C Reflection amplitude 

Here we discuss the variations of the Yang- Yang function with respect to the parameters U. 

Varying the action (13. 2p with the use of the on-shell condition ^Al^oi-fixed = 0, one easily 
derive the relation 

5 t A* = ttifa Sk , (C.l) 
where the constants fji can be regarded as regularized values of the solution fj at the boundary 

fji = lim [fi{w, w) — 2 li log \w — Wi\) . (C.2) 

\w— u>j\— >0 

It should be stressed that unlike li, which are "input" parameters applied with the problem, the 
values of the constants fji are not prescribed in advance but determined through the solution, 
i.e. it is rather part of the "output". 
Let us define 

6, = 2?f- 7 r*i)e*, (C.3) 
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where we use the notation 



2 2 V 2 

Bearing in mind the definition of the Yang- Yang function, (IC.lj) can be rewritten in the form 



Pi = m + I = £ (k + I ) • (C.4) 



^- = \o g e l . (c.5) 

As it follows from the small- p expansion (13.211) . 

&i = S(pi\p3+Pk)S(p i \p j -p k ) ex p($^^ p4n ) ' ( C - 6 ) 

n=l ^* 

Here k) stands for any permutation of (1, 2, 3) and 

, (r) , n) _ (p\ ~ - r(i +Pt + 9 )r(i +Pt -g) r(i-2 Pi ) r(i+jg) rr 
UJ r(|-p J + g)r(i-p 4 -g) r(i + 2 Pj ) r(i-f ) ' 

Notice that, with the parameter identifications ( 11.27ft . fll.28[t . 

2tt\ — 



(f ) ° l S( Pi \q) (CI 



R 

coincides with the so-called "reflection amplitude"- an important characteristic of the Fateev 
model (see Ref. [38J for details). 
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